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ABSTRACT 


Independent Earth orientation measurements taken by the space-geodetic techniques of 
lunar and satellite laser ranging, very long baseline interferometry, and the global positioning 
system have been combined using a Kalman filter The resulting combined Earth orientation 
series, SPACE2004, consists of values and uncertainties for Universal Time, polar motion, and 
their rates that span from September 28, 1976, to January 22, 2005, at daily intervals and is 
available in versions whose epochs are given at either midnight or noon. The space-geodetic 
measurements used to generate SPACE2004 have then been combined with optical astrometric 
measurements to form two additional combined Earth orientation series. (1) COMB2004, 
consisting of values and uncertainties for Universal Time, polar motion, and their rates that span 
from January 20, 1962, to January 22, 2005, at daily intervals and which is also available in versions 
whose epochs are given at either midnight or noon, and (2) POLE2004, consisting of values and 
uncertainties for polar motion and its rate that span from January 20, 1900, to January 20, 2005, at 
30.4375-day intervals. 
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INTRODUCTION 


Reference series of Earth orientation parameters (EOPs) obtained by combining 
independent measurements of the Earth’s orientation are generated annually at the Jet Propulsion 
Laboratory (JPL) in support of tracking and navigation of interplanetary spacecraft This report 
describes the generation of the most recent such combined Earth orientation series, SPACE2004, 
COMB2004, and POLE2004. Since the procedures used to generate these most recent series are 
similar to those used to generate previous such combinations, only a brief description of their 
generation is given here. Further details regarding the approach used at JPL to annually combine 
independent measurements of the Earth’s orientation can be found in Gross (1996, 2000) and 
Gross etal (1998). 


1 



SPACE2004 


Data Sets Combined to Form SPACE2004 

SPACE2004 is a combination of independent space-geodetic measurements of the Earth’s 
orientation. Table 1 lists the space-geodetic series used in generating SPACE2004, giving their 
identifiers, the number of measurements from each series that were actually incorporated into 
SPACE2004, and the time interval spanned by those measurements. Note that the University of 
Texas Center for Space Research (UTCSR) satellite laser ranging (SLR) Universal Time (UT) 
values were not used in generating SPACE2004 due to problems associated with separating this 
component of the Earth’s orientation from the effects of unmodeled forces acting on the satellite 
that cause the node of its orbit to drift (see Gross et al. 1998, p. 217 for further discussion about 
this point). For similar reasons, the Joint Center for Earth Systems Technology/Goddard Space 
Flight Center (JCET/GSFC) satellite laser ranging (SLR) length-of-day (LOD) values have not 
been used in generating SPACE2004 However, the JPL rapidly produced (quick-look) global 
positioning system (GPS) LOD series was used in generating SPACE2004 since (1) unlike the 
other available GPS LOD series, the JPL quick-look series is not cohstramed to UT1 
measurements and therefore consists of independent measurements of LOD, and (2) a stochastic 
model for the corrupting effects of the motion of the GPS satellite constellation on the JPL quick- 
look LOD values has been developed and is incorporated into the Kalman filter that was used here 
to combine the Earth orientation measurements. 

Since it was desirable to combine only independent measurements of the Earth’s 
orientation, only those series listed in Table 1 were used, even though other space-geodetic series 
are available from other analysis centers. When more than one series determined by the same 
measurement technique was used, care was taken to make sure that the measurements themselves 
were not included more than once In particular, polar motion measurements from the JPL GPS 
series were only used until the start of the International GPS Service (IGS) combined series 
EOP(IGS) 95 P 01 on January 1, 1995, polar motion measurements from the IGS combined series 
EOP(IGS) 95 P 01 were then used until the start of the IGS combined series EOP(IGS) 95 P 02 
on June 30, 1996; polar motion measurements from the IGS combined series EOP(IGS) 95 P 02 
were then used until the start of the accumulated IGS Solution Independent Exchange (SINEX) 
combined series EOP(IGS) 00 P 02 on February 27, 2000, and polar motion measurements from 
the accumulated IGS SINEX combined series EOP(IGS) 00 P 02 were used thereafter Similarly, 
measurements from the National Oceanic and Atmospheric Administration (NOAA) International 
Radio Interferometric Surveying (IRIS) Intensive UT1 series were used until it ended on December 
31, 1994, measurements from the USNO National Earth Orientation Service (NEOS) Intensive 
UT1 series were then used until it ended on December 4, 2000; and measurements from the GSFC 
NEOS Intensive UT1 series were used thereafter Finally, polar motion measurements from the 
UTCSR SLR series EOP(CSR) 96 L 01 were only used until the start of the JCET/GSFC SLR 
series on January 3, 1993, which was then used thereafter 
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Table 1. Data Sets Combined to Form SPACE2004 


Data Set Name Data Type 

Analysis Center Reference 

Data Span 

Number 

LLR (JPL05M01: VOL, UTO) 
McDonald Cluster 

LLR 

JPL 

Williams et al. (2002) 

Oct 5. 1976. to Dec. 31, 2004 

714 

OCA 

LLR 

JPL 

Williams et al. (2002) 

April 7 1984. to Jan. 19. 2005 

1048 

Haleakala 

LLR 

JPL 

Williams et al. (2002) 

Feb. 10. 1985, to Aug. 11 1990 

70 

UTCSR (CSR96L01) 
Lageos Polar Motion 

SLR 

UTCSR 

Eanes & Watkins (1996) 

Sep. 28, 1976, to Dec 30, 1992 

1846 

JCET/GSFC (JCET 04L83) 
Daily Lageos Polar Motion 

SLR 

JCET/GSFC 

Pavlis (2004) 

Jan. 3. 1993, to Jan. 2, 2005 

4365 

DSN (JPL97R01 T V) 
Califomia-Spain Cluster 

VLB I 

JPL 

Steppe et al. (1997) 

Nov 26, 1979. to Jan. 8, 2005 

836 

California- Australia Cluster 

VLBI 

JPL 

Steppe et al. (1997) 

Oct. 28,1978, to Sep. 30,1997 

698 

NOAA (NOAA95R02) 
IRIS Intensive UT1 

VLBI 

NOAA 

Ray et al. (1995) 

April 2, 1984. to Dec 31, 1994 

2394 

USNO (N9903) 

NEOS Intensive UT! 

VLBI 

USNO 

Eubanks et al. (1999) 

Jan. 4, 1995, to Dec. 4. 2000 

1497 

NASA/GSFC (GSFINT33) 
NEOS Intensive UT1 

VLBI 

GSFC 

NASA/GSFC VLBI Group (2005) 

Dec. 6, 2000. to Jan. 22, 2005 

904 

NASA/GSFC (GSF2005a) 
Multibaseline 

VLBI 

GSFC 

NASA/GSFC VLBI Group (2005) 

April 12. 1980, to Jan. 19, 2005 

3294 

NASA/GSFC (GSFC 1122) 
Westford-Fort Davis 

VLBI 

GSFC 

Ma & Gordon (1999) 

June 25, 1981, to Jan, 1 1984 

105 

Westford-Mojave 

VLBI 

GSFC 

Ma & Gordon (1999) 

March 21 1985, to Aug. 6, 1990 

18 

GPS (JPL 04APR05: LOD) 
JPL Rapid (Quick-Look) 

GPS 

JPL 

JPL Quick-Look Team (2005) 

Dec 27 1995, to Jan. 21. 2005 

3271 

GPS (JPL 21 APR04; Polar motion) 
Post-processed Flinn Analysis 

GPS 

JPL 

Heflin et al. (2004) 

June 10, 1992, to Dec 31 1994 

818 

GPS (IGS95P01: Polar motion) 
IGS Final Combined (ITRF93) 

GPS 

NRCan 

Kouba & Mireault (1997) 

Jan. 1 3 995, to June 29, 3 996 

546 

GPS (IGS95P02: Polar motion) 
IGS Final Combined (IGSOO) 

GPS 

CODE 

Mireault & Kouba (2000) 

June 30, 1996, to Feb. 26, 2000 

1334 

GPS (IGS00P02; Polar motion) 
IGS SINEX Combined (IGSOO) 

GPS 

NRCan 

Ferland (2004) 

Feb. 27 2000, to Jan. 22. 2005 

1789 


LLR, lunar laser ranging: JPL, Jet Propulsion Laboratory; VOL, variation of latitude; UT Universal Time: OCA Observatoire de la Cote 
d’Azur: UTCSR. University of Texas Center for Space Research: SLR, satellite laser ranging: JCL1 Joint Center for Earth Systems 
Technology: GSFC Goddard Space Flight Center: DSN, Deep Space Network: T transverse: V vertical; VLBI, very long baseline 
interferometry: NOAA, National Oceanic and Atmospheric Administration; IRIS, International Radio Interferometric Surveying, USNO. 
United States Naval Observatory: NEOS, National Earth Orientation Service, NASA, National Aeronautics and Space Administration; GPS, 
Global Positioning System; IGS, International GPS Service: ITRF International Terrestrial Reference Frame; NRCan, Natural Resources 
Canada; CODE, Center for Orbit Determination in Europe: SINEX, solution independent exchange. 

Data Preprocessing and Treatment of Tide-Induced Rotational Variations 

The Earth orientation series listed in Table 1 were first preprocessed by removing leap 
seconds from the UT1 values and, when necessary, by correcting the UT1 values to be consistent 
with the extended definition of Greenwich Sidereal Time (GST) as adopted by the International 
Earth Rotation Service (IERS; IERS 1997, p. 149). Since most of the series listed in Table 1 were 
already consistent with the extended definition of GST, this correction needed to be applied to only 
the NOAA IRIS Intensive UT1 senes. 
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Rotational variations caused by solid Earth and ocean tides were also removed from the 
UT1 values The effect of the solid Earth tides was removed by using the model of Yoder et al 
(1981), the model of Kantha et al. (1998) was used to remove the effect upon UT1 of the ocean 
tides at the Mf, Mf, and Mm tidal frequencies. Since the Yoder et al, (1981) model already 
includes a contribution from the equilibrium ocean tides, only the Kantha et al (1998) oceanic 
corrections to the Yoder et al. (1981) model were actually removed. Also note that the Kantha et 
al. (1998) model was used to remove the effect of ocean tides on only UT1, not on polar motion. 
Ocean-tide-induced polar motion variations were not removed from any of the polar motion 
observations, Finally, the only Earth orientation series listed in Table 1 that includes the effects of 
semidiurnal and diurnal ocean tides on the Earth’s orientation is the NOAA IRIS Intensive UT1 
series. This series included these effects by adding to the released UT1 values the model of 
Herring (1993, also see Herring and Dong 1994), Hence, the same Herring (1993) model was 
used to remove them. 

On June 30, 1996, the IGS reference frame changed from the International Terrestrial 
Reference Frame ITRF93 to ITRF94, on March 1, 1998, it changed from ITRF94 to ITRF96; on 
August 1, 1999, it changed from ITRF96 to ITRF97, on February 27, 2000, it changed from 
ITRF97 to the IGS realization of ITRF97 known as IGS97, on December 2, 200,1, it changed from 
IGS97 to the IGS realization of ITRF2000 known as IGS00; and on January 11, 2004, it changed 
from 1GS00 to a new IGS realization of ITRF2000 known as IGbOO, These changes in reference 
frames potentially introduce discontinuities into the IGS combined Earth orientation series. 
However, the IGS Final combined series EOP(IGS) 95 P 02 and the accumulated IGS SINEX 
combined series EOP(IGS) 00 P 02 used here have had each of these segments aligned to the same 
IGS reference frame. Thus, to within the uncertainty in determining the corrections required to 
align each segment, the IGS combined senes EOP(IGS) 95 P 02 and EOP(IGS) 00 P 02 used here 
should be reasonably homogeneous. They were, therefore, concatenated together with one common 
set of bias-rate corrections being determined for them, as described below 

Discontinuities can also be introduced into the JPL quick-look GPS LOD series by the 
above changes in the IGS reference frames or by other changes in processing strategy In fact, 
discontinuities at Modified Julian Date (MJD) 50956 (May 23, 1998), MJD 51391 (August 1, 
1999) and MJD 52383 (April 19, 2002) were empirically found in this series by differencing it with 
a reference combination of independent measurements. Therefore, the JPL quick-look GPS LOD 
series was partitioned into four segments with separate bias-rate corrections and a common 
uncertainty scale factor determined for each segment, as described below 

Adjustments Made to Space-Geodetic Series Prior to Combination 

Prior to combining the series listed in Table 1 to form SPACE2004, series-specific 
corrections were applied for bias and rate, the stated uncertainties of the measurements were 
adjusted by multiplying them by series-specific scale factors, and outlying data points were 
deleted. Values for the bias-rate corrections and uncertainty scale factors were determined by an 
iterative procedure wherein each series was compared, in turn, to a combination of all others, In 
order to minimize interpolation error (see Gross et al. 1998, pp 223-225), the comparison of each 
series to its reference combination was done at the epochs of the measurements of that series by 
generating its reference combination using a Kalman filter that interpolates to (and prints its EOP 
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estimates at) the exact epochs of those measurements. Also, both the bias-rate corrections and the 
uncertainty scale factors for all components of a given series were determined simultaneously in a 
multivariate approach using nonlinear weighted least squares. Using a multivariate approach 
allows the correlations between the components to be taken into account when determining the 
bias-rate corrections and uncertainty scale factors (see Gross et al. 1998, pp. 225). 

All the series listed in Table 1 were included in the iterative procedure except the JPL quick- 
look GPS LOD series, which was treated separately as described below Details of the iterative 
procedure, including: (I) the use of a reference series, SPACE2003 (Gross 2004), for initial bias- 
rate alignment, (2) the analysis of each data type in its natural reference frame; (3) the clustering of 
the McDonald LLR stations and, separately, the DSN VLBI stations in California, Australia, and 
Spain, (4) initial convergence on values for the series-specific bias-rate corrections and uncertainty 
scale factors prior to the start of outlier detection and deletion, and (5) final convergence on these 
values after detecting and deleting all data outliers, are described in Gross (1996, 2000) and Gross 
et al. (1998) and will not be repeated here. At the end of the iterative procedure, relative bias-rate 
corrections have been determined that make the series agree with each other in bias and rate; 
uncertainty scale factors have been determined that make the residual of each series have a reduced 
chi-square near one when differenced with a combination of all others, and outlying data points 
(those whose residual values are greater than four times their adjusted uncertainties) have been 
deleted. A total of 44 data points, or about 0 17% of those combined, were thus deleted. 

During the iterative procedure, the uncertainties of the JPL and IGS GPS polar motion 
series were adjusted in order to find and delete outlying data points, However, prior to combining 
these series with the others to form SPACE2004, their uncertainties were reset to their original 
values in order to avoid over-inflating them and hence overly smoothing the SPACE2004 series. 

The JPL quick-look GPS LOD series was not included in the above iterative procedure for 
bias-rate correction and uncertainty scale factor determination because of the use within the Kalman 
filter of a stochastic model for the uncertainties of this series. The stochastic model includes 
random walk and first-order autoregressive (AR1) components as well as a white noise component. 
These other components must be taken into account when adjusting the level of the white noise in 
the stochastic model. This was done after first separately determining and applying bias-rate 
corrections, but not uncertainty scale factors, to each of the four segments of the JPL GPS LOD 
series, These bias-rate corrections were determined by comparing each segment to a combination 
of all the other series after first applying to the other series the bias-rate corrections and uncertainty 
scale factors that were determined for them as described above. After applying the bias-rate 
corrections thus determined to each segment, an uncertainty scale factor and a bias-rate correction 
common to the entire JPL GPS LOD series was determined by iteratively comparing the JPL series 
to a combination of all of the series including the JPL series, In order for this reference 
combination to be quasi-independent of the JPL GPS LOD series, the Kalman filter that combined 
the series printed its estimate at the epoch of the LOD measurement before incorporating it into the 
state vector The uncertainty scale factor and bias-rate correction common to the entire JPL GPS 
LOD series were then determined that made its residual, when differenced with the reference series, 
have a reduced chi-square near one. The resulting uncertainty scale factor and bias-rate correction 
were then applied to the entire JPL GPS LOD series, and the process was repeated until the scale 
factor changed by less than 0.2%. Besides determining the uncertainty scale factor and bias-rate 
corrections, this process also led to the detection and deletion of 21 outlying data points. 
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Table 2. Adjustments to Space-Geodetic Series 


Data Set Name 

Bias, mas 


Rate, mas/yr Uncertainty Scale Factor 

LLR (JPL05M01) 

VOL 

UT0 


VOL UT0 

VOL 


UT0 

McDonald Cluster 

0.086 ± 0.090 

0.178 ±0.099 

-0.015 ±0.021 -0.014 ± 0.022 

1.471 


1372 

OCA 

-0.083 ± 0.048 

0.007 ± 0.038 

0.016 ± 0.007 -0.001 ±0.005 

2.399 


1.994 

Haleakala 

-1.498 ±11 19 

-1.438 ±0.745 

-0.275 ±0.218 -0.248 ±0.156 

2.652 


2.502 

DSN (JPL97R01) 

T 

V 


T V 

T 


V 

Califomia-Spain Cluster 

-0.739 ± 0.025 

-0.229 ± 0.057 

0.1 86 ±0.003 0.166 ±0.009 

1.644 


1.282 

California- Australia Cluster 

0.708 ±0.018 

0.502 ± 0.053 

0.041 ±0.008 -0.027 ±0.024 

1.488 


1.166 

NASA GSFC (1122) 

T 

V 


T V 

T 


V 

Westford-Fort Davis 

8.619 ± 3.591 

3.110 ±5.648 

0.790 ±0.362 0.237 ±0.566 

1.390 


0.886 

Westford-Mojave 

0.869 ±1.042 

0.099 ± 1 752 


0007 -0004 

2.368 


0.760 

GSFC (GSF2005a) PMX 

PMY 

UT1 

PMX 

PMY UT1 

PMX 

PMY 

UT1 

Multibaseline 0.257 ± 0.007 

0.109 ±0.006 

0.071 ±0.011 

-0.021 ±0.001 

0.012 ± 0.001 0.010 ±0.002 

1.825 

1.654 

1.756 

NOAA (95R02) PMX 

PMY 

UT1 

PMX 

PMY UT1 

PMX 

PMY 

UT1 

IRIS Intensive 


0.280 ± 0.023 


0012 ±0.007 



1.044 

USNO (N9903) PMX 

PMY 

UT1 

PMX 

PMY UT1 

PMX 

PMY 

UT1 

NEOS Intensive 


0.009 ±0.041 


0.014 ±0008 



1.323 

GSFC (GSFINT13) PMX 

PMY 

UT1 

PMX 

PMY UT1 

PMX 

PMY 

UT1 

NEOS Intensive 


0.044 ± 0.300 


0.009 ±00.10 



1.270 

GPS (JPL Rapid LOD) PMX 

PMY 

LOD 

PMX 

PMY LOD 

PMX 

PMY 

LOD 

MJD 50078-50956 


0.302 ±0.148 


0.031 ±0.034 



1.614 

MJD 50956-51391 


-0.840 ± 0.455 


0090 ±0.076 



1.614 

MJD 51391-52383 


-0.495 ± 0.224 


0.111 ±0028 



1.614 

MJD 52383-53391 


0.433 ± 0.223 


0.022 ±0.021 



1.634 

UTCSR (96L03) PMX 

PMY 

LOD 

PMX 

PMY LOD 

PMX 

PMY 

LOD 

Lageos -0.21 1 ± 0.030 

0.771 ±0.024 


0097 ± 0.008 

0.100 ± 0007 

1.060 

0.864 


JCET (JCET 04L83) PMX 

PMY 

LOD 

PMX 

PMY LOD 

PMX 

PMY 

LOD 

Daily Lageos -0.353 ± 0.006 

0.822 ± 0.006 


-0.001 ± 0.001 

0.094 ± 0.001 • • 

2.497 

2.208 


GPS (JPL 21APR04) PMX 

PMY 

LOD 

PMX 

PMY LOD 

PMX 

PMY 

LOD 

JPL Post-Flinn -0.009 ± 0.0 1 4 

-0.0 M ±0.013 


-0005 ± 0.014 

0.001 ±0.012 

3.000 

1.000 


GPS (IGS95P01) PMX 

PMY 

LOD 

PMX 

PMY LOD 

PMX 

PMY 

LOD 

Final Combined -0.024 ± 0.051 

-0.004 ± 0.045 


0.008 ± 0018 

0.001 ±0016 

1.000 

1.000 


GPS (IGS95P02) PMX 

PMY 

LOD 

PMX 

PMY LOD 

PMX 

PMY 

LOD 

Final Combined 0 .020 ± 0.0 1 0 

0.017 ± 0.008 


0001 ±0.001 

0.002 ±0.001 

1.000 

1.000 


GPS (IGSO0P02) PMX 

PMY 

LOD 

PMX 

PMY LOD 

PMX 

PMY 

LOD 

SINEX Comb. 0.020 ± 0.0 1 0 

0.017 ±0.008 


0001 ±0.001 

■0.002 ±0003 

LOGO 

3.000 



Reference date for bias-rate adjustment is 1993.0. See Table 1 footnotes, mas, milliarcseconds; PMX. polar motion X; PMY polar motion Y 


Finally, each series was placed within a particular IERS reference frame by applying to it 
an additional bias-rate correction that is common to all the series. This additional correction was 
determined by first combining all the series after applying to them the relative bias-rate corrections 
and uncertainty scale factors that had been determined for them as described above. This 
intermediate combination was then compared to the IERS combined Earth orientation series 
EOP(IERS) C 04 (IERS 2004, Section 3.5 1) during the interval 1987-2004 in order to obtain the 
additional bias-rate correction required to make it, and therefore each individual series, agree in bias 
and rate with the EERS series. 

The total bias-rate corrections and uncertainty scale factors determined by the procedures 
outlined above are given in Table 2. The values of the bias-rate corrections given in Table 2 are the 
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sum of: (1) all the incremental corrections applied during the iterative procedure, (2) the corrections 
applied to initially align the series with each other, and (3) the additional, common correction 
applied in order to place each series within that particular IERS reference frame defined by the 
IERS combined Earth orientation series EOP(IERS) C 04 The values of the uncertainty scale 
factors given in Table 2 are the products of all the incremental scale factors determined during the 
iterative procedure. The uncertainties of the bias-rate corrections given in Table 2 are the la 
standard errors of the incremental bias-rate corrections determined during the last iteration. There 
are no bias -rate entries in Table 2 for components that were either not used (e.g., the UTCSR SLR 
UT1 component) or not available (e.g., the NOAA IRIS Intensive polar motion components), 
Note that the same IERS rate correction is applied to all the data sets, including those such as the 
Westford-Mojave single baseline VLBI series for which no relative rate correction could be 
determined. Therefore, the rate corrections given in Table 2 for those series for which no relative 
rate corrections could be determined are simply the IERS rate correction, but given in the natural 
reference frame of that series. In these cases, uncertainties for the rate corrections are not given. 
Also note that the entries for the bias-rate corrections m Table 2 for the IGS combined series 
EOP(IGS) 95 P 02 and EOP(IGS) 00 P 02 are the same. Since these series were initially given 
within the same IGS reference frame, they were merged together and a common bias-rate correction 
was determined for them. The entries in Table 2 for the uncertainty scale factors of the JPL and all 
three IGS GPS polar motion series are unity since, as discussed above, the uncertainties of these 
series have been reset to their original values prior to forming SPACE2004, 


Combined EOP Series: SPACE2004 

A Kalman filter was used to combine the series listed in Table 1 after the bias-rate 
corrections and uncertainty scale factors listed in Table 2 had been applied to them, The resulting 
combined Earth orientation series, SPACE2004, consists of values (Figure 1) and la standard 
errors (Figure 2) for polar motion. Universal Time, and their rates spanning September 28, 1976, to 
January 22, 2005, at daily intervals and is available in versions whose epochs are given at either 
midnight or noon. Leap seconds have been restored to UT1, and the effects of the solid Earth and 
ocean tides have been added back to the UT1 values using the same models for these effects that 
were originally used to remove them from the raw series: namely, the Yoder et al. (1981) model for 
the solid Earth tides and the Kantha et al. (1998) model for the ocean tides. However, semidiurnal 
and diurnal ocean tidal terms have not been added to and are therefore not included in the 
SPACE2004 UT1 values. 

Figure 3 is a plot of the difference between the SPACE2004 polar motion, UT1, and LOD 
values and those of the IERS combined Earth orientation series EOP(IERS) C 04. These two 
series are very consistent with each other, especially after March 1, 1998 when the root-mean- 
square (rms) of their difference is only 0.065 milliarcseconds (mas) for the x-component of polar 
motion, 0.060 mas for the y-component, 0.031 milliseconds (ms) for UT1, and 0.022 ms for 
LOD Prior to 1984, the difference between these two series exhibits greater variability and even 
some systematic behavior, particularly in the x-component of polar motion. This systematic 
behavior is due to differences in the approaches used here and by the IERS to correct the bias and 
rate of the individual series before they are combined. 
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Figure 2. Plots of the l a formal errors of the x-component of polar motion (2a), y-component of polar motion (2b), UT1-UTC (2c), and excess length-of-day 
(2d) as given by the combined Earth orientation series SPACE2004, The insert within panels 2a, 2b, and 2c shows that component's post- 1984 uncertainties on 
an expanded scale with the same units: milliarcseconds (mas) for polar motion, milliseconds (ms) for UT1 UTC. 
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Figure 3 .Plots of the difference between the IERS combined Earth orientation series EOP(IERS) C 04 and SPACE2004 formed by subtracting the SPACE2004 
values from those of the IERS series, The difference between the x-component of polar motion is shown in 3a, the difference between the y-component is 
shown in 3b, the difference between UT1-UTC is shown in 3c, and the difference between the excess length-of-day is shown in 3d, 







COMB2004 


COMB2004 extends SPACE2004 by additionally incorporating the optical astrometric 
polar motion and UT1 series that was determined at the Bureau International de l’Heure (BIH) 
from an analysis of time and latitude observations by Li (1985, also see Li and Feissel 1986), This 
BIH optical astrometric series consists of values and uncertainties for polar motion and UT1 that 
span from January 5.0, 1962 to December 3 1 .0, 1981 at 5-day intervals. 

Data Preprocessing and Treatment of Tide-Induced Rotational Variations 

The BIH optica] astrometric series was first preprocessed by removing leap seconds from 
the UT1 values and by correcting the UT1 values to be consistent with the extended definition of 
GST, as adopted by the IERS (IERS 1997, p. 149), Rotational variations caused by solid Earth 
and ocean tides were also removed from the UT1 values. The same models that were used to 
remove the tidal effects from the series combined to form SPACE2004 were also used to remove 
them from the BIH series, namely, the Yoder et al, (1981) model for the solid Earth tides and 
the Kantha et al (1998) model for the Mf, Mf, and Mm ocean tides, However, since the BIH UT1 
measurements represent an average value over a 5-day-long observation window, and since 5 days is 
a substantial fraction of the monthly and shorter-period tides, the amplitudes of these solid Earth and 
ocean tidal terms were attenuated prior to their removal from the BIH UTl measurements. (See 
Gross 1996, p. 8735 and Gross et al. 1998, pp. 226-227 for further discussion about this point.) 


Adjustments Made to BIH Series Prior to Combination 


The preprocessed BIH optical astrometric series was combined with the space-geodetic 
series that comprise SPACE2004 after first (1) correcting the BIH series to have the same bias, 
rate, annual, and semiannual terms as SPACE2004, (2) applying a constant multiplicative scale 
factor to the measurement uncertainties of the BIH series so that its residual, when differenced with 
SPACE2004, had a reduced chi- square of one; and (3) deleting those data points, if any, whose 
residual values were greater than four times their adjusted uncertainties. Due to software 
limitations associated with the need to correct the annual and semiannual terms of the BIH series, 
the above adjustments were determined separately for each component of the BIH series in a 
univariate approach rather than simultaneously in a multivariate approach as was done for the series 
combined to form SPACE2004. The procedure used to determine these bias-rate and seasonal 
term corrections, uncertainty scale factors, and outlying data points has been described before 
(Gross 1996, pp. 8735-8738) and will not be repeated here. The annual and semiannual terms of 
the BIH series were adjusted in order to correct for systematic, seasonally varying effects that are 
known to be present in optical astrometric measurements, Since the values of both the BIH and 
SPACE2004 series are given at midnight, interpolation error (see Gross et al 1998, pp. 223-225) is 
automatically minimized when differencing these two series for the purpose of determining the 
adjustments to be made to the BIH series. Tables 3 and 4 give the resulting uncertainty scale factors 
and values and lo standard errors of the corrections to the bias, rate, annual, and semiannual 
terms thus determined for the BIH series. When determining these uncertainty scale factors and 
the corrections to the bias, rate, and seasonal terms, no outlying data points were detected. 
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Table 3. Adjustments to Bias, Rate, and Stated Uncertainty of Optical Astrometric Series 


Data Set 

Bias, mas 



Rate, mas/yr 


Uncertainty Scale Factor 


PMX 

PMY 

UT1 

PMX 

PMY 

UT1 

PMX 

PMY 

UT1 

BIH 

-1.878 ±0.856 

0.527 ±0.716 

-8.214 ±0.770 

0.889 ±0.481 

1.308 ± 0.192 

5.386 ±0.319 

1.807 

1.682 

1.928 

ILS 

-48.346 ± 2.384 

-3.210 ± 1.688 


0.261 ±0.487 

0.052 ± 0.344 


2.185 

1.522 



Reference date for bias-rate adjustment of BIH series is 1980.0. Reference date for bias-rate adjustment of ILS series is 1970.0. 


Table 4. Adjustment to Annual and Semiannual Terms of Optical Astrometric Series 


Data Set Coefficient of Sine Term, mas Coefficient of Cosine Term, mas 


PMX PMY UTJ PMX PMY UT1 


BIH 

annual 

semiannual 


-6.352 ± 1.023 
2.058 ± 1.043 


-5.408 ± 0.695 
0.076 ±0.698 


5.154 ±0.813 
-0.272 ± 0.828 


-2.916 ± 1.087 9.860 ±0.787 -0.912 ± 0.873 

0.983 ± 1,062 0.597 ± 0.732 1.667 ± 0.856 


ILS 

annual 

semiannual 


-3.051 ±3.327 
0.030 ± 3.335 


10.067 ±2.347 
9.356 ±2.354 


8.002 ± 3.363 -10.206 ± 2.376 

2.217 ±3.351 1.390 ±2.367 


Reference date for adjustment of annual and semiannual terms of BIH series is 1980.0. 
Reference date for adjustment of annual and semiannual terms of ILS series is 1970.0, 


Combined EOP Series: COMB2004 

A Kalman filter was used to combine the BIH senes with the adjusted space-geodetic series 
that comprise SPACE2004 after first applying to the BIH series the uncertainty scale factors and 
corrections to the bias, rate, annual, and semiannual terms that are given in tables 3 and 4. The 
resulting combined Earth orientation series, COMB2004, consists of values (Figure 4) and 1 
standard errors (Figure 5) for polar motion, Universal Time, and their rates that span from January 
20, 1962, to January 22, 2005, at daily intervals and is available in versions whose epochs are given 
at either midnight or noon Leap seconds have been restored to UT1 and the effects of the solid 
Earth and ocean tides have been added back to the UT1 values using the same models for these 
effects that were originally used to remove them, namely, the Yoder et al (1981) model for the solid 
Earth tides and the Kantha et al. (1998) model for the ocean tides. The full amplitude (i.e., no tidal 
terms attenuated) of the effects of the solid Earth and ocean tides at the epoch of the time tag were 
added back to the UT1 values. Semidiurnal and diurnal ocean tidal terms have not been added to 
and are therefore not included in the COMB 2004 UT1 values, 
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A COMBINED EARTH ORIENTATION SERIES COMB2004 
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Figure 4. Plots of the x-component of polar motion (4a), y-component of polar motion (4b), UT1-UTC (4c). and excess length-of-day (4d) as given by the 
combined Earth orientation series, COMB2004 The discontinuous changes in the plot of UT1-UTC are caused by the presence of leap seconds. Prior to the 
introduction of leap seconds in 1972, Coordinated Universal Time (UTC) was adjusted by introducing step and rate changes designed to keep it close to UT1 
(e.g., IERS 1997, Table II-3), the effect of which is also readily apparent in 4c Note that the UT1 -UTC values displayed in 4c include the tidal variations, 
whereas the excess length-of-day values shown in 4d do not 








x-component of polar 
series COMB2004, 







POLE2004 


No optica] astrometric observations taken at the stations of the International Latitude 
Service (ILS) were used when creating the BIH optical astrometric series used in COMB2004 (Li 
1985, Li and Feissel 1986), The ILS polar motion measurements (Yurrn and Yokoyama 1980), 
which are based solely upon latitude observations made at the ILS stations are, therefore, 
independent of those comprising COMB2004 and have therefore been combined with them to form 
POLE2004 Being based solely upon latitude observations, the ILS series contains no UT1 
measurements, but consists solely of polar motion measurements that span 1899 8-1979.0 at 
monthly intervals Although no uncertainties are given with the individual polar motion values, the 
precision with which they have been determined is estimated to be 10-20 mas (Y umi and Yokoyama 
1980, p. 27), An initial uncertainty of 15 mas was, therefore, assigned to each of the ILS polar 
motion values, Since this assigned measurement uncertainty will be adjusted later, its initial value is 
arbitrary as long as it is not zero and serves merely as an a priori estimate to be used in the series 
adjustment procedure described below 

The ILS series was combined with COMB 2004 to form POLE2004 after - (I) correcting the 
ILS series to have the same bias, rate, annual, and semiannual terms as COMB2004, (2) applying a 
constant multiplicative scale factor to the measurement uncertainties of the ILS series so that its 
residual, when differenced with COMB 2004. had a reduced chi -square of one, and (3) deleting 
those data points, if any, whose residual values were greater than four times their adjusted 
uncertainties. Again, due to software limitations associated with the need to correct the annual and 
semiannual terms, these adjustments were determined separately for the x- and y-components of the 
ILS polar motion series by fitting a bias, rate, and these seasonal terms to the difference of the ILS 
series with COMB2004 during 1962.0 to 1979.0, The measurement uncertainties of the ILS polar 
motion values were adjusted by determining and applying a scale factor that made the residual of 
this fit have a reduced chi-square of one. During this procedure to determine uncertainty scale 
factors and bias, rate, and seasonal term corrections, no outlying ILS data points were deleted since 
no data points had residual values greater than four times their adjusted, uncertainties. Tables 3 and 
4 give the resulting uncertainty scale factors and values and la standard errors of the corrections to 
the bias, rate, annual, and semiannual terms thus determined for the ILS series. 

A Kalman filter was then used to combine the ILS series with the adjusted BIH and space- 
geodetic series that comprise COMB 2004, after applying to the ILS series the uncertainty scale 
factors and corrections to the bias, rate, annual, and semiannual terms that are given in Tables 3 and 
4. The resulting combined Earth orientation series, POLE2004, consists of values (Figure 6a and 
6b) and la standard errors (Figure 6c and 6d) for polar motion and its rate that span from 
January 20, 1900, to January 20, 2005, at 30.4375-day intervals. 
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DISCUSSION 


The Kalman filter that was used here to combine Earth orientation measurements contains a 
stochastic model of the process that is used to propagate the state vector and its associated state 
covariance matrix forward in time to the epoch of the next measurement. For polar motion 
excitation, the stochastic model includes a random walk term with equal noise forcing both the x- 
and y-components of excitation and having a white noise power spectral density of 246.4 mas 2 /day 
(Morabito et al 1988), This level of polar motion excitation process noise m the Kalman filter was 
recently increased to 739.2 mas 2 /day in order to better match the observed spectrum of polar 
motion excitation. Increasing the excitation process noise reduces the level of smoothing applied to 
the polar motion components of the propagated state vector and increases the covariance associated 
with those components. Thus, the SPACE2004, CQMB2004, and POLE2004 polar motion and 
polar motion excitation values are not as heavily smoothed as are those of previous combinations 
produced at JPL, such as those of SPACE2003, and the uncertainties assigned to the SPACE2004, 
COMB2004, and POLE2004 polar motion and polar motion excitation values are somewhat larger 

Since a Kalman filter has been used in generating SPACE2004, C.OMB2004, and 
POLE2004, the resulting polar motion and UT1 values are smoothed to a degree that depends upon 
both the spacing between the measurements being combined, which determines how far the state 
vector and state covariance matrix must be propagated forward in time, and the uncertainties that 
have been assigned to the measurements. Since improvements to the observing systems, both in 
the hardware and software and in the number of systems, have led to increasingly precise 
determinations of the Earth’s orientation, and since the time resolution of the measurements has 
generally increased in concert with the measurement precision, the degree of smoothing applied to 
the SPACE2004, COMB2004, and POLE20G4 values is a function of time, with the earlier values 
being more heavily smoothed than the more recent values, 

Daily EOP values are reported in SPACE2004 since the NOAA IRIS, USNO NEOS, and 
GSFC NEOS Intensive UT1 values are given at daily intervals, as are the GPS and JCET/GSFC 
SLR values (although gaps exist in each of these data sets). However, prior to the start of these data 
sets, the measurements combined to form SPACE2004 are given less frequently, therefore, the 
Kalman filter used to combine these measurements also interpolates them in order to produce a 
series of equally spaced values. In order to be consistent with SPACE2004, daily EOP values are 
also reported in COMB2004 even though the BIH optical astrometric series used in COMB2004 is 
given at 5-day intervals. Thus, SPACE2004, COMB 2004, and POLE2004 are equally spaced 
series of smoothed, interpolated Earth orientation parameters. 

The combined Earth orientation series SPACE2004, COMB2004, and POLE2004 are 
available from: (1) JPL’s Space Geodetic Science and Applications Group by anonymous ftp to 
<ftp://eulerjpl.nasa.gov/keof/combinations/2004>, (2) NASA’s Crustal Dynamics Data Informa- 
tion System (CDDIS) by anonymous ftp to <ftp://cddis.gsfc.nasa.gov/products/jpl/2004>, or (3) 
the author upon request to Richard.Gross@jpl.nasa.gov 
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